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Abstract 

The state of the searches for isolated black holes, non-accreting black holes in 
binary systems and, finally, the accreting black holes in the X-ray binaries is pre- 
sented. The third category is, by far, the most important source of information 
about the Galactic black holes and it is the main topic of this review. Different 
proposed signatures of an accreting black hole are reviewed (including the grow- 
ing evidence for the presence of the event horizon in some systems). The list of 
51 black hole candidates compiled with the help of these signatures is presented. 
To obtain a conclusive evidence that the candidate is really a black hole, one 
needs a dynamical mass estimate indicating that the mass of the compact object 
is greater than ~ 3Mq. The observational methods of obtaining such estimates 
are briefly discussed. The application of these methods leads to the list of 18 
objects that could be considered confirmed black holes. Their estimated masses 
are found to fall in the range 3.5 -^IGM©. Finally, some general remarks concern- 
ing relation of black hole candidates to different classes of X-ray binaries such 
as X-ray Novae, microquasars and high/low mass X-ray binaries are presented. 
The implications of the properties of black hole binaries for the theory of stellar 
evolution are also discussed. 



1 Introduction 

Half a century ago, a black hole was just an exotic theoretical concept. Thirty years 
ago, the first Galactic black hole candidate (Cyg X— 1) was proposed (Bolton, 1972). 



For the next ten years, Cyg X— 1 remained the only candidate and some researchers 
put a lot of effort to devise a model of this system without invoking a black hole. With 
the advent of the subsequent black holes in binary systems (LMC X— 3, Mon X— 1), 
the motivation for such efforts substantially diminished. By today, black holes became 
a well established and extensively investigated class of observed astrophysical objects 
and their existence is no longer seriously disputed. Recently, searches for isolated 
Galactic black holes started and at least one object (microlens MACHO-99-BLG-22 = 
OGLE-1999-BUL-32) became a strong candidate. 

In the initial part of this article, I shall briefly describe the state of the searches for 
isolated black holes and for non-accreting black holes in binary systems. Most of my 
review will be devoted to black hole candidates in the X-ray binaries. I shall describe 
the methods of detecting black holes in such systems, review the properties of the 
candidates (51 objects) and discuss briefly some implications of these properties for 
the stellar evolution theory. 

2 Isolated Black Holes 

The evidence for detection of isolated black holes is much weaker than the one for binary 
black holes. However, some candidates were proposed among the objects responsible 
for long microlensing events and also among the objects listed in the Sloan Digital Sky 
Survey catalogues. It seems, that only one reasonably convincing case for an isolated 
black hole was demonstrated so far. 

2.1 Extended Microlensing Events 

The long time scale of a microlensing event may be attributed either to a large mass 
of the lens or to the small relative transverse velocity of the lens with respect to the 
lensed object. The long events show, usually, the magnification fiuctuations, reflecting 
the motion of the Earth. This effect permits to calculate the so called "microlensing 
parallax" which is a measure of the relative transverse motion of the lens with respect 
to the observer. Assuming standard model of the Galactic velocity distribution, one 
is then able to perform a likelihood analysis, which permits to estimate the distance 
and the mass of the lens. Basing on such analysis, several long events detected by the 
MACHO and OGLE projects were suggested as, possibly, caused by black hole lenses 
(Bennett et al. 2001, Mao et al. 2001, Bennett et al. 2002, Smith et al. 2002a). Agol 
et al. (2002) indicated that such analysis might overestimate the number of black hole 
lenses. At the moment, only one of the proposed candidates (MACHO-99-BLG-22 = 
OGLE-1999-BUL-32), with the time scale ~ 600 days) remains as a robust candidate 
for a black hole lens. 

2.2 Isolated Accreting Black Holes 

It has been estimated (Agol and Kamionkowski, 2001) that substantial number of iso- 
lated black holes may accrete from the interstellar medium with the efficiency sufficient 
to be detectable in deep X-ray surveys. Models of the emission from such objects sug- 
gest that they should also produce a synchrotron radiation detectable in the visible 



wavelengths (Ipser and Price 1982, McDovell 1985). Chisholm et al. (2002) searched 
through the 3.7 million objects in the Sloan Digital Sky Survey data and selected 150 
000 objects with colors and properties consistent with synchrotron visible spectrum. 47 
of these objects are ROSAT X-ray sources. For seven of these sources optical spectra 
excluding the candidacy (2 quasars and 5 F type stars) were observed. The remaining 
40 sources are listed by the authors as possible black hole candidates. 

3 Non- Accreting Black Holes in Binary Systems 

Discovery of a black hole in a binary system with no mass transfer (and so no substan- 
tial X-ray emission) is rather difficult. Some W-R stars have unseen companions of 
several solar masses. However, these companions do not have to be compact objects. 
They might be just normal main sequence stars, unseen because their luminosities are 
much lower than the luminosities of the bright W-R primaries. Cherepashchuk (1998) 
presents the list of 17 single-line W-R binaries, suspected of containing unseen compact 
components. Some of these components may be black holes, but no strong conclusions 
can be made. 

4 Signatures of an Accreting Black Hole in a Binary 
System 

Over the years, the researchers have accumulated a number of signatures, indicating 
that an accreting compact object in an X-ray binary is a likely black hole. These 
signatures could be ordered into several major classes: 

• Signatures in the X-ray spectrum 

• Signatures in the X-ray variability 

• Level of the X-ray luminosity 

• Evidence for the presence of the event horizon 

4.1 Signatures in the X-Ray Spectrum 

Over the last quarter of a century, a number of characteristics were proposed as X- 
ray signatures of an accreting black hole. Some of them, like flickering or bimodal 
behaviour (state transitions) were found also in neutron star systems, so they are 
no longer beheved to be unique properties of black holes. Those, listed below, should 
be considered as useful indicators (but not a definite proof) of the presence of a black 
hole. 

(1) Presence of an ultrasoft component (kT~ 0.1—0.3 keV) in the X-ray 
spectrum. This component, dominating in the high^ (soft) state, is believed to 
originate from the inner region of an accretion disc. This region has larger radius 
and lower temperature in the case of a black hole than in the case of a neutron star 
(see the point (4) below) and therefore produces softer X-ray spectrum. 

^"high" and "low" refer to X-ray flux below 10 keV, and not necessarily to the luminosity integrated 
over the entire X-ray range 



(2) Presence of a power law component (extending to several tens or several 
hundreds keV) in the X-ray spectrum. This component, presumably produced 
through the up-comptonization of soft photons and sometimes described as an 
ultrahcird component, dominates in the low (hard) state. It is present also in 
some neutron star systems, but then it is not as hard as in black hole systems. In 
hard states, black holes are harder than neutron stars (Sunyaev, 2001). 

(3) Absence of the characteristic ~ 2 keV blackbody component at all levels 
of X-ray emission (and, in particular, at high X-ray luminosity: Lx ^ 10^^ erg 
s~^). This component, prominent in the spectra of non- pulsating neutron star 
sources at high luminosity levels, is believed to originate on the surface of a weakly 
magnetized neutron star. 

Both the absence of the ~ 2 keV blackbody component and the presence of the supersoft 
component (see the point (1) above) contribute to the fact that black holes have softer 
X-ray colors than the neutron stars and, in the soft states, both classes of objects 
occupy different regions in the X-ray color-color diagram (Done and Gierlihski, 2002). 
They occupy different regions in the color-color diagram also in hard states, since then 
the power law component is harder in black holes than in neutron stars (point (2) 
above) . 

(4) Lcirge inner radius of the disc and low color temperature of the inner 
edge of the disc. These parameters arc estimated from multi-color blackbody fits 
to the ultrasoft component of the X-ray spectrum. The results of such fits indicate, 
that the inner radii of discs around black hole candidates (BHCs) are, consistently, 
by a factor ~ 4 larger than for the discs around neutron stars (NSs) (Tanaka and 
Lewin, 1995). They indicate also, that the inner parts of discs around BHCs are 
cooler than the discs around NSs. 

4.2 Signatures in the X-Ray Variability 

Relatively recently, some characteristics of X-ray variability were proposed as signatures 
of an accreting black hole. They include: 

(1) Relatively narrow power density spectrum with high frequency cut-off at 
~ 10 Hz, during hard/low states. Sunyaev and Revnitsev (2000) noticed that 
there are substantial differences between hard/low states power spectra of BHC 
and NS systems: (1) NSs have generally significantly broader spectra and (2) they 
demonstrate substantial power at high frequencies up to ^ 500 Hz, while BHCs 
have little power above v ~ 10-50 Hz. 

(2) Low frequency (0.1—3 Hz) quasi-periodic oscillations (QPOs). Both BHC 
and NS systems show low frequency QPOs. Generally, the low frequency QPOs of 
BHCs have lower frequencies (0.1—3 Hz), than those of NSs (5—60 Hz) (see e.g. 
Psaltis and Norman, 2000). However, it should be noted, that some BHCs show 
QPOs with frequencies ~ 5—10 Hz (and so within NSs range). 

(3) High frequency quasi-periodic oscillations. High frequency QPOs (or kHz 
QPOs) are seen both in NS systems (200 - 1330 Hz) and in BHC systems (41 
- 450 Hz). There are different proposed mechanisms to explain their origin. At 
present, the leading contenders are: the parametric cpicyclical resonance in the 
inner accretion disc (Abramowicz and Kluzniak, 2001) or the different modes of 
oscillations of the inner disc (e.g Wagoner et al., 2001). The fact that in some 



systems (and especially in some BHCs) the observed QPO frequencies form integral 
harmonics and that the frequencies seem to scale with the mass like (see 
Remillard et al., 2002 for the most recent review) gives strong support to the first 
solution. In both theories, the values of the frequencies of the pair QPOs provide 
us with some constraints on the mass and the angular momentum of the compact 
object. Let us note, that (similarly as for low frequency QPOs) the frequencies 
are generally lower for BHCs than for NSs (presumably because of scaling with the 
mass - see Abramowicz et al., 2000). High frequency QPOs are therefore a good 
diagnostic tool to distinguish between accreting BHs and accreting NSs. 

4.3 Level of the X-Ray Luminosity 

A new class of point X-ray sources, called "ultraluminous compact X-ray sources" 
(ULXs) was estabhshed in recent years (Makishima et al., 2000). Objects of this type 
reside in many nearby galaxies (although none was found in our Galaxy or in M31). 
ULXs have properties of typical X-ray binaries, except that their X-ray luminosities 
are much higher: Lx ~ 10^^ -r 10^° erg s~^). This is too much for an accreting neutron 
star (unless it is a super- Eddington source). The most straightforward interpretation 
is that ULXs are accreting black holes with the masses ~ 10 -^ lOOOM©. However, the 
alternative explanation, assuming strongly coUimated emission, is also possible and, so 
far, there is no consensus on this point. 

4.4 Evidence for the presence of the event horizon 

(1) The lack of Type I bursts from some X-ray Novae. X-ray Novae (XNe, called 
also Soft X-Ray Transients or SXTs) change their accretion rates (as inferred from 
the X-ray luminosities) by many orders of magnitude. Therefore, there must be 
time intervals, when the accretion rates fall into certain critical range of values, 
corresponding to the presence of the Type I bursts (thermonuclear flashes due to 
unstable nuclear burning of the accreted matter on the surface of the NS). About 
50 XNe are known at present. In roughly 85 % of them, the compact object is a 
BHC, in the remaining ~ 15 % it is a NS. In most of the NS XNe, Type I bursts 
are observed from time to time. However, not a single Type I burst was ever seen 
in any of the BHC XNe. This is a strong argument in favour of the presence of 
the event horizon, since Type I bursts can occur only if the solid surface is present. 
Narayan and Heyl (2002) investigated the stability of the nuclear burning on the 
surfaces of accreting compact objects of 1.4 and 10 Mq for different accretion rates 
(assuming the presence of the solid surface in each case). In agreement with many 
earlier investigations, they found instability for some accretion rates in the case of 
1.4 Mq object (corresponding to the NS). But they found instability also for the 10 
Mq object. It seems, therefore, that independently on the mass. Type I bursts must 
appear, if only solid surface is present. The lack of these bursts implies, therefore, 
the presence of the event horizon around the compact object. 

(2) Very low quiescent luminosities of some X-ray Novae. Between the out- 
bursts, XNe retain some (although very low) level of X-ray luminosity, which in- 
dicates that some marginal accretion is still present. This, so called "quiescent 
luminosity" is so low, that it was not detectable, until the recent generation of the 



X-ray satellites. The comparison of the quiescent luminosities of BHC XNe and 
NS XNe shows a very substantial difference between the two classes of objects. For 
comparable orbital periods (which implies comparable sizes of the orbits and, so, 
the comparable accretion rates), the BHC XNe are, systematically, by about two 
- two and a half orders of magnitude dimmer than NS XNe (Narayan et al. 2001, 
Garcia et al. 2001). The most obvious explanation of this difference is the presence 
of the solid surface (where most of the accretion energy is released) in the NS XNe. 
BHC XNe are "blacker" because they have event horizons. 

5 Black Hole Candidates in X-Ray Binaries 

Table 1 contains the list of BHCs in the X-ray binaries, compiled with the help of 
criteria discussed in the previous paragraph. I removed Cyg X-3 from the list, as there 
are huge (two orders of magnitude!) discrepancies in the estimates of its mass function 
(Schmutz et al., 1996, Hanson et al., 2000). However, I believe that Cyg X-3 will return 
to the hst in not too distant future. The hst contains, at present, 51 objects. Thanks 
to RXTE and Bcppo-SAX it grew recently at a rate of 2—3 new BHCs per year. 

I should stress that the signatures used to compile the list serve only as useful diag- 
nostic tools: they help us to select systems suspected of harboring a black hole. They 
constitute, therefore, very useful criteria for compiling a list of black hole candidates. 
To obtain a conclusive evidence (and so, to include an object into the list of confirmed 
black holes), we have to obtain a dynamical mass estimate of the compact object and 
to demonstrate that this mass is too high for a neutron star. 

6 Confirmed Black Holes (Candidates with Dynam- 
ical Mass Estimates) 

6.1 Upper Mass Limit for Neutron Stars 

• Theory 

The maximum permitted mass of a neutron star (so called Oppenheimer— Volkoff 
mass) depends on the equation of state of matter at the densities higher than the 
nuclear density. Unfortunately, so far, this equation is not well known. The theory is 
not able to discriminate between different proposed models. For seemingly reasonable 
equations of state, the value of Oppenheimer— Volkoff mass Mqv ranges from ~ 1.5Mq 
for soft equations of state (density grows relatively fast with pressure) to ~ 2.7 Mq 
for hard equations of state (density grows relatively slowly with pressure). If one 
assumes only general relativity as the correct theory of gravitation and the principle of 
causality (the speed of sound not greater than the speed of light) then, for arbitrarily 
exotic equation of state (constrained only by the requirement of the causality), one 
obtains Mqv ~ 3.2Mq for non-rotating configuration and Mqv ~ 3.9M0 for maximally 
rotating configuration. 

• Observations 

The most precise mass determinations are available for neutron stars that are mem- 
bers of binary radio pulsar systems composed of two neutron stars. The radio pulsar 



B1913+16 and its NS companion have measured masses 1.4411 and 1.3874 Mq, re- 
spectively, with the error ± 0.00035 Mq. For the system B1534-I-12 the masses of both 
components are 1.339 Mq with the error ± 0.003 Mq. These two systems set the range 
of the observed values of the neutron star masses. There are numerous other determi- 
nations: three other binary NS systems, nine systems composed of radio pulsar and 
white dwarf or main sequence star and eleven accreting X-ray pulsars (van Paradijs 
1998, Thorsett and Chakrabarty 1999). The results range from 1.06 to 1.88 Mq, but 
the errors are much larger, so that there is no determination that would be inconsistent 
with the range estabhshed from the two systems mentioned above. The general rule is 
such that the more precise is the determination, the closer it approaches the value ~ 

1.4 Mq. 

To summarize: all observational determinations of neutron star masses are consis- 
tent with the values in the range 1.34 1.44M0. 

• Conclusion 

If the accreting component of an X-ray binary has mass greater than ~ 3M0, it 
must be a black hole. 

6.2 Dynamical Mass Estimates 

There are several observational parameters used to estimate the mass of the compact 
component. They arc briefly discussed below. 

• The mass function 

This is the most important observational parameter used to constrain the mass 
of the compact component. The mass function f{Mx) is calculated from the radial 
velocity variations of the optical companion: 



where Kopt is the semiamplitude of the radial velocities of the absorption lines of the 
optical component (in km/sec) and P is the orbital period (in days). The mass function 
is related to the masses of both components by: 



where M^ and Mopt arc the masses of the compact and the optical components and 
i is the inclination of the orbit. The value of f{Mx) gives an absolute lower limit to 
the compact component mass. In many cases, this is enough to conflrm the presence 
of a BH in the system (if /(Mj.) is ^ 3Mq). Out of 18 BHs with the mass estimates, 
7 systems have mass functions ^ 5Mq and 3 other systems have f{M^) in the range 
3 4Mq. However, in some cases, we need more precise constraints. To obtain them, 
we have to estimate the mass ratio q = Mopt/Mx and the inclination of the orbit i. 

• The rotational broadening of the absorption lines of the optical com- 
ponent 

Prom the measurements of the lines, we determine the projected rotational velocity 
at the equator of the optical component Vj-otsini. Assuming the corotation of the 
optical component with the orbital motion (for a Roche lobe filling component, it is a 
very good assumption), we have: 





/(M,) = MlsinH/{M^t + M,f 



(2) 



Vrotstm = OAQKopt{q/{l + g))'/' (3) 

With the help of this equation, one can determine the mass ratio q. The intrinsic 
width of the absorption hnes is very small ( ~ few km/sec), while the typical rotational 
broadening is of the order of few tens km/sec and can be, relatively easily, measured. 
This method has been applied, with a substantial success, to many BHC systems (see 
Orosz, 2002). 

• The radial velocities of the emission lines of the accretion disc 

The shifts of these lines reflect the orbital motion of the compact component and 
so permit us to determine the mass ratio directly: 

q = Kem/Kopt (4) 

where Kem is the semiamplitude of the radial velocities of the emission lines. This 
method of determining the mass ratio is completely independent from the previous 
one. Unfortunately, it is rather uneasy to implement, since the emission lines are very 
broad (~ 2000 km/sec), while the orbital shifts are of the order of few tens km/sec. 

• The amplitude of the ellipsoidal light Vciriations 

Due to filling of the Roche lobe, the optical component is tidally distorted and due 
to rotation it exhibits the ellipsoidal light variations (double sinusoid per orbital period 
if the rotation is synchronous) . The amplitudes of these variations in V and I are given 
by: 

A\/ = 0.26smV(l + g) (5) 

A/ = 0.24smV(l + g) (6) 
As may be seen, the dependence on the mass ratio q is rather weak (for most of the 
BH XNe q falls in the range 0.05 -r- 0.2) and, therefore, the ellipsoidal light variations 
provide us with a valuable information about the inclination i (even, if we do not know 
q). In practice, the procedure is not that simple, since the optical light is usually (even 
in the quiescence) contaminated by the residual contribution from the accretion disc. 
Using light curves in different colors, one obtains different values for the inclination 
(see the discussion of A 0620—00 in section 6.3). To make the situation worse, this 
residual contribution is frequently variable. Using the observational data from different 
time intervals, one again obtains different values for the inclination. Extraction of the 
true value requires, often, a very careful modeling (Froning and Robinson 2001, Gelino 
et al. 2001b). 

• The mass-spectral type relation for the optical component 

The typical optical components of BHXNe are the lower main sequence stars, which 
satisfy reasonably well the mass-spectral type relation. This relation may be used to 
estimate the mass of the compact component, if the mass ratio is unknown. Since the 
mass of the optical component is usually quite small (below 1 Mq), even substantial 
uncertainty (say factor of two) does not influence significantly the mass estimate for 
the compact component. 

• High frequency QPOs 

As discussed in section 4.2 p. (3), the sets of frequencies of the kHz QPOs provide 
us with some information on the mass of the compact object. 



6.3 Black Holes Candidates with Dynamical Mass Estimates 



The dynamical mass estimates, indicating the mass of the compact component greater 
than ~ SMq, are now available for 18 X-ray binaries. The estimated masses of black 
holes in these 18 systems are given in Table 1. These estimates are based on the liter- 
ature and discussed briefly below. Recently, Orosz (2002) performed new calculations 
of all masses using the published values of mass functions, inclinations and rotational 
velocities and employing a simple Monte Carlo code. Such approach has an advan- 
tage of a relatively uniform treatment of all systems and a disadvantage of omitting 
important pieces of information considered by different authors while discussing the 
individual systems. 

Below, I shall briefly discuss the parameters of these 18 systems. 
(1) Cyg X-1 

At present, there are no doubts about the presence of black hole in the system. 
However, in spite of three decades of investigations, there is still substantial uncertainty 
concerning the masses of both components. The value of the mass function given by 
Gies and Bolton (1982): /(M^) = 0.25 ± O.OIM© was only shghtly modifled recently 
by Brocksopp et al. (1999): /(M^) = 0.244 ± O.OOSM©. Gies and Bolton (1986a, 
1986b), analyzing the emission lines of the stellar wind, the rotational broadening of 
the absorption lines of the optical component and the photometric V-band light curve, 
found that this component must be close to flUing its Roche lobe (the flUing factor must 
be greater than 0.9). Prom the same analysis, they estimated the inclination of the 
orbit to be i = 33° ± 5°. From polarimetric measurements (in three colors) Dolan and 
Tapia (1989) found a rather wide range of values for inclinations: i = 62°(+5°, —37°). 
Paczyhski (1974) has shown that basing only on • mass function • lack of X-ray eclipses 
• obvious fact that the optical component cannot be larger than its Roche lobe and • 
the spectral type and the photometry of the optical component, one can obtain lower 
limits to the masses of both components as functions of the distance to the system. 
Unfortunately, this distance is not well established. However, most likely, Cyg X— 1 
belongs to Cyg OB— 3 association and so its distance is d = 2.15 ± 0.07 kpc (Massey 
et al., 1995). With this distance, the present day cahbrations of effective temperatures 
and bolometric corrections: — 32000 K (Herrero et al., 1995), B.C. = -3.18 (Vacca 
et al., 1996) and Massey's photometry, following Paczyhski's reasoning, one obtains 
~ 8Mq and ~ 28M0 as lower limits to the masses of the compact and the optical 
components, respectively (Ziolkowski, 2002). 

Quite independent constraints on the masses may be obtained from analysis of the 
evolutionary status of the optical component (Ziolkowski, 2002). At the distance of 
2.15 ± 0.07 kpc, its luminosity must be Mboi ~ —9.4 ± 0.2. It can be argued that it 
must be a core hydrogen burning configuration and then its mass cannot be smaller 
than ~ 31±3M0. This lower limit is consistent with the values discussed above. Even 
if the distance to the system was as small as 1.8 kpc (the smallest value quoted in 
the literature), then still the optical component should be more massive than ~ 24 
Mq. Altogether, the evolutionary analysis indicates that the mass of HDE 226868 is 
probably in the range 30 -^ 40 Mq and the initial mass of its progenitor was probably 
45 -7- 55 Mq. The mass of the compact component is probably in the range 15 17 
Mq. These results values are fully consistent with the values M^pt = 33 ± 9Mq and 
Mx = 16 ± 5Mq, obtained by Gies and Bolton (1986a) from the essentially distance 



independent (but model dependent) analysis. 
Orosz (2002) obtained = 6.85 H- 13.25M0. 

(2) LMC X-3 

The mass function is given by Cowley et al. (1983): /(M^) = 2.3 ± O.3M0. The 
same authors estimated the rotational broadening of the absorption lines of the optical 
component to be 130 ± 20 km/sec and concluded that the inclination of the orbit must 
be limited to i ^ 50° -j- 70° (the upper limit comes from the lack of X-ray eclipses). 
Their estimate of was 7 IAMq with the most likely value ~ 9Mq. Kuiper et al. 
(1988) included the photometric data (exhibiting the eUipsoidal hght variations) and 
obtained i = 64° -i- 70° and = 4.5 ^ 6.5M0. 

Orosz (2002), using rotational velocity of Cowley et al. and inclination of Kuiper 
et al., obtained = 5.94 9.17Mq. 

(3) LMC X-1 

The mass function, given by Hutchings et al. (1987), is f(M^) = 0.14 ± 0.05Mo. 
The lack of X-ray echpses implies i ^ 60°. Hutchings et al. obtain ~ 4 -^ IOMq 
for the mass of the compact component. Shrader and Titarchuk (1999) argue for 
much higher value {M^ ^ I6M0), basing on their modeling of the accretion disc X-ray 
spectrum. It is not clear how to reconcile these two estimates. Orosz (2002) quotes 
the value of Hutchings et al. 

(4) SS 433 

This celebrated object has been extensively reviewed on many occasions (see e.g. 
Margon, 1984). Here, I shall consider only the data important for the mass estimates. 

The system displays two basic periods: the binary orbital period (Porb ~ 13.087 ± 
0.003 d) and the period of the precession of the accretion disc and jets (Pprec ~ 162.5 
d). The inclination of the orbital plane (from the kinematic model of the jets) is 
i = 78.82° ± 0.11°, the precession angle 6 = 19.80° ± 0.18°, the distance to the object 
(from the radio observations of the jets) d ~ 5 kpc. The visual extinction is estimated to 
be Ay ~ 8 magnitudes. The system is very luminous in the optical band: Lopt > 2x10^^ 
erg s~^. About 80 % of this luminosity comes from the photosphere of geometrically 
thick accretion disc, which can be treated as a flattened star. The observed He II 
emission lines also originate on the surface of the disc and reflect its orbital motion. 
Until recently, no spectral lines produced by the second member of the binary (the 
mass donor or the "optical component") were seen. Therefore, untypically, only the 
mass function /(Mopt) could be estimated. Unfortunately, its value is not estabhshed 
very well. Looking through the literature, one can flnd /(Mopt) ~ IO.IMq (Crampton 
and Hutchings, 1981), /(Mopt) = 7.7(+3.0, -2.4)Mq (Fabrika and Bychkova, 1990) 
and /(Mopt) = 1-90 ± O.25M0 (D'Odorico et al., 1991). Fabrika and Bychkova gave 
the most convincing arguments supporting their value and, therefore, it will be used 
in further considerations. Substantial efforts were undertaken to derive the mass ratio 
Q — Mx/Mopt from the modeling of the optical light curves. The light curves in two 
colors are two-dimensional (B and V magnitudes are modulated with two periods: Porb 
and Pprec) and rather complicated. The modeling assumes mutual partial eclipses of 
two photospheres (that of the accretion disc and of a "normal star", which must be a 
massive early type star) with possible addition of a hot spot on the surface of the disc. 
Two such analysis carried out by Leibowitz (1984) and Antokhina and Cherepashchuk 
(1985), under different assumptions about the geometrical shape of the disc gave very 
similar results. Both papers got thickness of the disc equal approximately to its radius 



(ratio of polar to equatorial radii of about 0.5). Lcibowitz obtained Q > 0.8 and the 
preferable value of Q ~ 1.2. Antokhina and Cherepashchuk got acceptable solutions 
for Q 0.4 and Q <^ 1.1 (they preferred Q ~ 1.2). Quantitatively, the results of these 
modeling were very unprecise (with the mass function of Fabrika and Bychkova, the 
value of Afx is in the range 4.2 62Mq), but they confirmed the presence of a black 
hole in the system. 

Quite recently, Gies et al. (2002) observed, finally, the absorption lines of the 
optical component (mass donor). They estimated the spectrum as A-type supergiant 
and measured the semiamplitude of the radial velocities: K^pt — 100 ± 15 km/sec 
(which corresponds to f{Mx) = 1.36(+0.71, — O.52)M0. Taking the semiamplitude of 
the radial velocities of the emission lines from Fabrika and Bychkova: = 175 ± 20 
km/sec, they estimate the mass ratio Mx/Mopt = 0.57 ± 0.11 and the masses Mopt = 
19 ± 7Mq and = 11 ± 5Mq. 

(5) GRO J0422+32 

The most recent values of the mass function, are given by Harlaftis et al. (1999): 
/(M^) = 1.13±O.O9M0 and Webb etal. (2000): /(M^) = 1.191±O.O21M0. Filippenko 
et al. (1995) estimated the mass ratio from the radial velocities of the accretion disc 
emission hues and got q = 0.109±0.009. This value is consistent with the ones obtained 
from the rotational broadening of the absorption lines of the optical component: q — 
0.116(+0.079, -0.071) (Harlaftis and Filippenko, 2000) and q = 0.111(+0.048, -0.022) 
(Webb ct al., 2000). Taking into account f{Mx) and q and assuming normal (for its 
spectral type) mass of the secondary (~ O.4M0) Filippenko et al. (1995) obtained 
i = 48° ± 3° and = 3.57 ± O.34M0. 

Orosz (2002) obtained = 3.66 ^4.97M0. 

(6) A 0620-00 

The mass function is /(M^) = 2.72±O.O6M0 (Marsh et al., 1994). The same authors 
used the absorption lines of the optical component to estimate the mass ratio: q = 
0.067 ±0.010. Mass ratio was estimated also from the radial velocities of the accretion 
disc emission lines: Haswell and Shaffer (1990) obtained q = 0.094(+0.022, -0.015), 
while Orosz et al. (1994) got q = 0.074 ± 0.006. The determinations of the inclination 
of the orbit, estimated from the ellipsoidal light variations have a long and complicated 
history. Haswell et al. (1993), using variations in the UBVR bands, found i = 70°±3.5°, 
which lead to the mass of the compact component — 3.8 4.4M0. Shahbaz et al. 
(1994a) used the variations in the JK bands and got i — 37° ± 5°, which lead to 
Mx = 14 ± 7M0. The second determination seemed to be closer to the truth, since 
the flux in the infrared bands should be less contaminated by the residual contribution 
of the accretion disc. However, the next determination, using the H band photometry 
(Froning and Robinson, 2001), did not solve the problem. The authors found, first, that 
the light curves change from year to year and, second, that the determination of i is 
strongly model dependent and admits a wide range of values for the inclination. They 
decided only to conclude that i should be in the range 38° ^ 75°, which corresponds 
to Mj; = 3.3 I3.6M0. Finally, Gelino et al. (2001b) succeeded in modeling their 
new JHK photometry (they had to introduce cold spots on the surface of the optical 
component) and found i = 40.8° ± 3° and = 11.0 ± 1.9M0. 

Orosz (2002) obtained = 8.70 ^ I2.86M0. 
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NOTES: 

Porfe — orbital period 

Opt. Sp — optical spectrum 

X-R — X-ray variability 

C — comments 

Mbh— mass of black hole component 
Ref — references 

H optical counterpart was identified, but the spectrum was not obtained 

IR — optical coTintcrpart was seen only in infrared 

(O) — name of the optical object 

T — transient 

RT — recurrent transient 

pers — persistent 

/xQ — microquasar 
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(7) GRS 1009-45 

The mass function, based on observations with the help of Keck II telescope, is 
/(Ma,) = 3.17 ± O.12M0 (Filippenko et al., 1999). The same observations yielded 
precise measurements of the radial velocities of the accretion disc emission lines, which 
permitted to calculate the mass ratio q = 0.137 ± 0.015. The lack of X-ray eclipses 
implies i ^ 80°. Assuming that the secondary has a normal mass for its spectral type 
(~ O.6M0), the authors obtain ~ 4.4Mq and i ~ 78°. The mass of the compact 
object cannot be smaller because of the limits set by the lack of X-ray eclipses. It can 
be somewhat larger (up to ~ 4.7M0 if the mass of the secondary is ~ O.75M0). 

Orosz (2002) quotes, after Filippenko et al. (1999), = 3.64 ^ AMMq. 

(8) XTE J1118+480 

The most recent value of mass function, given by Wagner et al. (2001) is /(Ma,) = 
6.1 ± O.3M0. McClintock et al. (2001) found ellipsoidal hght variations in the I band 
with the amplitude ~ 0.15 mag. Assuming that the secondary has a normal mass 
for its spectral type (~ 0.5Mq), and considering different possible values of the disc 
contribution to the I-band emission (40% to 66%), McClintock ct al. found i = 52° -=-80° 
and Mx = 7.2 ^ 13.2M0. Wagner et al. (2001) analyzed their R band photometry and 
found ellipsoidal hght variations with the amphtude ~ 0.2 mag. Modehng their hght 
curves and permitting variable contribution from the disc (decreasing with time from 
72% to 64%), they found i = 81° ± 2° and = 6.0 ^ 7.7M0 for plausible optical 
component masses in the range 0.09 -r O.5M0. 

Orosz (2002) obtained = 6.48 H- 7.19M0. 

(9) GS 1124-684 

The mass function, given by Orosz et al. (1996), is /(M^) = 3.01 ± O.15M0. In 
an earlier paper, Orosz et al. (1994) found (from the rotational broadening) q = 
0.133 ± 0.019. Using these data and including the BVRI photometry (ellipsoidal light 
variations) the authors obtained i = 54° 65° and = 5.0 -r 7.5Mq. 

Later estimates were somewhat more conservative but they did not differ substan- 
tially: f{M^) = 3.34±O.15M0, q = 0.128±0.04 (Casares et al., 1997), and i = 39°^74°, 

= 3.8 ^ IO.5M0 (Shahbaz et al, 1997). 

The most recent analysis (Gelino et al., 2001a) uses the new JK photometry. They 
get i = 54° ± 2° and = 6.95 ± O.6M0. 

Orosz (2002) obtained = 6.47 ^ 8.I8M0. 



(10) 4U 1543-475 

The first mass function, given by Orosz et al. (1998), was f{M^) = 0.22±0.02Mq. 
Orosz et al. analysed also ellipsoidal light variations seen in the V and I hght curves. 
Assuming that the secondary has a normal mass for its spectral type (~ 2.3 2.6Mq), 
they found i = 30° ± 2° and M^^5±1Mq. 

In the most recent paper, Orosz (2002) quotes for the mass function f{M,j.) = 0.25± 
O.OIM0, for the rotational velocity i = 46 ± 2 km/sec and for the inclination 

i = 20.7° ± 1.5°. He concludes that = 8.45 ^ 1O.39M0. 

(11) XTE J1550-564 

The mass function, given by Orosz et al. (2002) is f{M^) = 6.86 ± O.71M0. 
The same authors obtained the tentative measurement of the rotational broadening 
{vrotsini = 90 ± 10 km/sec) and determined (from the ellipsoidal light variations) the 
orbital inclination: i — 73° ± 2°. Prom careful modeling of all observational data, they 
get Mx = 9.68^11.58M0. 

Orosz (2002) quotes another solution from Orosz et al. (2002): Mx = 8.36 
1O.76M0. 

(12) GRO J1655-40 

The most recent mass function of this famous microquasar was given by Shahbaz 
et al. (1999): /(M^) = 2.73±O.O9M0. Israelian et al. (1999) measured the rotational 
broadening: Vrotsini = 93 ± 3 km/sec), which implied the mass ratio q = 0.40 ± 0.03. 
From the careful modeling of the ellipsoidal light variations obtained with the BVIJK 
photometry, Greene et al. (2001) get i = 70.2° ±1.2° and = 6.3±O.3M0. Shahbaz 
(2002) from precise modeling of the spectrum of the optical component got q = 0.419 ± 
0.028 and = 5.99 ± O.42M0. 

Orosz (2002), practically, repeats the solution of Greene et al. (M^ = 6.03 
6.57M0). 

(13) H 1705-250 

The mass function, given by Harlaftis et al. (1997) is f{M^) = 4.65 ±O.21M0. The 
same authors got q ~ 0.014(+0.019, —0.012) (from the rotational broadening of the 
optical component lines). Filippenko et al. (1997) obtained q ~ 0.01 from the emission 
lines of the disc, but they did not consider this result to be realistic. Analysis of the 
ellipsoidal light variations in B and V bands leads to i = 70° ± 10° (Remillard et al., 
1996), while variations in R band give i ^ 48° ^ 51° (Martin et al., 1995). 

Remillard et al, (1996) obtained = 6 ± IM0, while Filippenko et al. (1997) 
got Mj, = 6.4 6.9M0 (assuming the mass of the optical component in the range 
0.3 H- 0.6Mq). Harlaftis et al. (1997) state carefully that = 4.9 H- 7.9M0, but their 
analysis seems to support rather the upper part of this range (5.7 -^ 7.9M0). 

Orosz (2002) obtained 714 = 5.64 ^ 8.3OM0. 

(14) XTE J1819-254 

The mass function, given by Orosz et al. (2001), is f{M^) = 2.74 ± O.12Af0. The 
same authors measured the rotational broadening and found v^-otsini = 124 ±4 km/sec. 
They analyzed also the ellipsoidal light variations and got i = 60° -j- 71°. Finally, they 
concluded that the mass of the compact component is in the range 8.73 11.7OM0. 

In the most recent paper, Orosz (2002) quotes for the mass function /{M^) — 
3.13 ± O.I3M0, for the rotational velocity Vrotsini = 98.9 ± 1.5 km/sec and for the 
inchnation i = 75° ± 2°. He concludes that = 6.82 ^ 7.42Mq. 



(15) XTE J1859+226 

The preliminary mass function was determined by Filippenko and Chornock (2001): 
f{Mx) ~ 7.4 ± 1.1M0. Ellipsoidal light variations in R band were analyzed by Zurita 
et al. (2002), who found only that the inchnation must be high. The rotational 
broadening was not measured yet. Assuming that the mass of the optical component 
is ^ IM©, one finds ^ 8.0 ^ 1O.2M0. 

Orosz (2002) obtained = 7.6 12.0Mq. 

(16) GRS 1915+105 

The mass function of the celebrated microquasar was determined by Greiner et al. 
(2001): /(Mj.) = 9.5 ± 3.0Mq. The authors assumed then that the inclination is given 
by the orientation of the jets: i = 70° ± 2°. Assuming, in turn, that the mass of the 
optical component is ^ I.2M0, they obtained = 14 ± AMq. The same value is 
quoted by Orosz (2002). 

(17) GS 2000+251 

The mass function, based on observations with the help of Keck I telescope, is 
f{Mx) = 5.01 ± 0.12Mq (Harlaftis et al., 1996). The same observations yielded the 
estimate of the rotational broadening {vrotsini = 86 ± 8 km/sec), which permitted to 
determine the mass ratio q — 0.042 ± 0.012. Analysis by Callanan et al. (1996), who 
included photometry (ellipsoidal light variations) in JK bands and assumed normal 
(for its spectral type) mass of the secondary (~ 0.4 O.9M0), gave i — 65° ± 9° 
and Mx = 8.5 ± 1.5Mq. The more recent analysis, incorporating the new J-band 
photometry (Leeber et al., 1999) gives i ~ 75° and ~ 6.55M0. 

Orosz (2002) quotes i = 64.0°±1.5° for the inclination and gets = 7.15^7.78M0. 

(18) GS 2023+338 

The mass function, given by Casares and Charles (1994), is f{M,j,) = 6.O8±O.O6M0. 
The same authors analyse the rotational broadening ( = 39.1 ± 1.2 km/sec), 

obtaining for the mass ratio the value q = 0.060 ± 0.005. Shahbaz et al. (1994b) 
included photometry (eUipsoidal light variations) in K band and found i — 56° ± 4° 
and Mx ~ I2M0 (10 I5M0). Sanwal et al. (1996) analyzed photometry in H band 
and found that cannot be larger than ~ 12.5M0, and most likely ^ 9.5M0. 

Orosz (2002) quotes i — 38.8° ± 1.1° for the rotational broadening and gets Mx = 
10.06 H-13.38M0. 

7 General Remarks 

• Brief Statistics 

The list of BHCs (see Table 1) contains 51 objects (including 18 confirmed BHs). 
About 44 of them (including 13 confirmed) have low mass companions (and so are 
members of low mass X-ray binaries), while only about 7 (5 confirmed) arc members of 
high mass X-ray binaries. Let us recall, that for the X-ray pulsars (strongly magnetized 
neutron stars) the distribution is just reversed: out of more than 80 accreting X-ray 
pulsars, presently known, only few are found in low mass X-ray binaries. 

• Black holes and X-ray Novae 

Galactic binary black holes show clear preference for a specific class of X-ray bi- 
naries, namely, for X-ray Novae. This is true also vice versa: great majority of X-ray 
Novae contain black holes. Out of ~ 49 XNe known so far (the precise number depends 



on the definition of an X-ray Nova, so perhaps it is better to talk about SXTs), only 
in ~ 7 the compact object is a neutron star. The remaining systems (~ 42) contain 
BHCs. Since typical XNe remain dormant most of the time, the true number of binary 
black holes might be easily by one or two orders of magnitude larger than the, 51 
objects, listed in Table 1. Most of these black holes reside in presently dormant XNe. 
Let us recall, that XNe containing BHCs are now being discovered at a rate of 2—3 
new objects per year. 

• Black holes and microquaseirs 

Out of 16 known galactic microquasars, two contain neutron stars ( Sco X— 1 and 
Cir X— 1), in one the nature of the compact object is unclear (Cyg X— 3) and in 13 the 
compact object is a BHC Among these 13, there are three superluminal jet sources. 
There is some evidence that two of superluminal sources contain Kerr (rapidly rotating) 
BHs. 

• Black holes and stellair evolution 

Analysis of parameters of BHC binaries and their possible evolutionary histories 
(Kalogera, 1999a,b; Wellstein and Langer, 1999; Brown et al., 1999; Nelemans et al., 
1999, Belczyhski and Bulik, 2002) leads to interesting conclusions concerning stellar 
evolution. Among the most important are: (1) BHs are formed more easily than we 
used to believe; the limiting main sequence mass Mns/bh above which the star ends 
as a black hole is ~ 20-^25 Mq rather than ~ 40—50 M©, as claimed previously; (2) 
during the supernova explosion of a helium star less than one half of the stellar mass 
is ejected, (3) the mass loss due to stellar winds from helium stars might be smaller 
than predicted by canonical formulae. 
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